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Abstract 

Objectives: To confirm genotype diversities of clinical isolates of Bordetella 
pertussis and to evaluate the risk of pertussis outbreak in Korea. 
Methods: Seven housekeeping genes and 10 antigenic determinant genes from 
clinical B. pertussis isolates were analyzed by Multilocus sequence typing (MLST). 
Results: More variant pattern was observed in antigenic determinant genes. 
Especially, PtxS1 gene was the most variant gene; five genotypes were observed 
from eight global genotypes. In the bacterial type, the number of observed 
sequence types in the isolates was seven and the most frequent form was type 1 
(79.6%). This major sequence type also showed a time-dependent transition 
pattern. Older isolates (1968 and 1975) showed type 1 and 6 in housekeeping 
genes and antigenic determinant genes, respectively. However, these were 
changed to type 2 and 1 in isolates 1999-2008. This transition was mainly 
attributed to genotype change of PtxS1 and Fim3 gene; the tendency of genotype 
change was to avoid vaccine-derived genotype. In addition, there was second 
transition in 2009. In this period, only the sequence type of antigenic determi- 
nant genes was changed to type 2. Based Upon Related Sequence Types (BURST) 
analysis confirmed that there were two clonal complexes (ACCI and ACQ I) in the 
Korean isolates. Moreover, the recently increased sequence type was revealed as 
AST2 derived from AST 3 in ACCI. 

Conclusions: Genotype changes in Korean distributing strains are still progressing 
and there was a specific driving force in antigenic determinant genes. Therefore 
continuous surveillance of genotype change of the distributing strains should 
be performed to confirm interrelationship of genotype change with vaccine 
immunity. 
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1. Introduction 

Pertussis is an infectious respiratory disease that is 
especially prevalent in newborn babies and infants who 
have not completed the basic vaccine schedule [1]. The 
pathogen of pertussis has been known as Bordetella 
pertussis since 1906 [2]; it is disseminated to only 
humans by infectious droplets released from carriers 
[3]. In the initial stage of infection, symptoms akin to 
a common cold appear, followed by characteristic 
symptoms of pertussis, such as paroxysmal cough and 
whoop sounds in breathing. After this paroxysmal stage, 
patients usually recover, with a long convalescent 
stage accompanied by coughing [1,4]. However, some 
patients may develop more severe symptoms such as 
pneumonia, convulsion, encephalopathy, rib fracture, 
and even death [4]. 

The global incidence rate of pertussis has greatly 
decreased following the introduction of a vaccine. 
However, there are several considerations in control of 
pertussis incidence. One is that local outbreaks occur 
even in highly vaccinated countries [5]; another is the 
transition of the infected age group from infants to adults 
[6]. One of the inferred reasons is that distributing strains 
evade the protection effect of vaccination by alteration 
of their surface antigens [7—9]. Currently, genetic 
polymorphisms were reported in four major antigenic 
determinant genes, such as pertussis toxin, filamentous 
hemagglutinin, pertactin, and fimbriae [10—15], which 
encode protein components of acellular pertussis 
vaccine. It was also reported that the distributing strains 
show distinct genotype from vaccine strain [9,16]. A 
third inferred reason is waning of vaccine-induced 
immunity with increasing age [1,9,17]. It is known that 
the vaccine's efficacy starts to decrease after 6—8 years 
from the last vaccination. Therefore, vaccine-induced 
immunity is maintained for approximately 12—15 
years, then immunity gradually diminishes [1,18]. 

In Korea, the average pertussis incidence rate from 
2000 has been maintained at 11.5 cases per year [19]. 
However, in 2009, the number of pertussis patients 
rapidly increased to 66 cases, a 500% increment. Most 
patients (85.2%) were infants aged <1 year [20]. 
Although the number of patients decreased to 27 in 
2010, it was still higher than the annual average for 
2000—2008. Furthermore, 36.6% of contacted persons 
from eight cases were confirmed as PCR-positive, and 
most (80%>) were adults. 

Considering the above, analysis of the reason for 
this transient increase of pertussis incidences in 2009 
was clearly needed, so as better to control pertussis 
outbreaks. Therefore we confirmed the genotypes in 
antigenic determinant genes and housekeeping genes of 
the clinical B. pertussis strains. Then, the genotype 
variations were analyzed by multilocus sequence typing 
(MLST). 



MLST is a sequence-based typing method to identify 
isolates into clusters defined as sequence type (ST) with 
an identical genotype set [21]. The major advantages 
of the MLST method in bacterial typing are to allow 
rapid generation of clear results and easy comparisons 
between STs. Usually, seven core metabolic house- 
keeping genes originating from multilocus enzyme 
electrophoresis (MLEE) analysis are used in MLST 
analysis. However, other genes apart from housekeeping 
genes are also applicable if there are low genetic vari- 
ations in housekeeping genes [22]. 

In this report, we analyzed genotype variations of 
seven housekeeping genes and 1 0 antigenic determinant 
genes to evaluate the genotype variation level of domes- 
tically distributed strains. In addition, the tendency of 
genotype change was analyzed to confirm the current 
status of pertussis incidence and to produce reliable data 
for pertussis control. 

2. Materials and Methods 

2.1. Bacterial isolates and culture condition 

A total of 74 B. pertussis strains, mainly isolated 
between 2000 and 2009 in Korea, were analyzed 
(Figure 1). All isolated strains were preserved in deep 
freezer (— 80°C) until analysis. Beside these clinical 
isolates, additionally five clinical strains isolated before 
2000 and 3 reference strains (ATCC 10380, ATCC 
9797, Tohama I) were included in all experiments. 
Reagan Lowe medium containing 1 0% horse blood was 
used to culture the isolated strains; inoculated medium 
was incubated in a humid incubator at 37°C for 5—7 
days. 

2.2. Isolation of genomic DNA 

Genomic DNA of the isolates was prepared by 
QIAamp DNA mini kit (QIAGEN) according to the 
manufacturer' s instructions . 

2.3. PCR of target genes and sequencing 

PCR primers and PCR condition of most test genes 
were cited from previous reports [12,23—25] and the 
Bordetella Multi Locus Sequence Typing website 
(http://pubmlst.org/bordetella/) (Table 1). In addition, 
new primers for FHA, OmpQ, BapC, CyaA, Vag8, and 
TcfA genes were applied. After confirming amplification 
of target genes by specific primer sets, the PCR products 
were purified by QIAquick PCR purification kit (QIA- 
GEN) and their nucleotide sequences confirmed by 
direct sequencing method. 

2.4. Determination of genotype profiles and 
completion of sequence types 

To determine the genotypes of housekeeping genes, 
sequence information of tested strains was queried to 
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5 
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KNT65 


KNT51 
KNT53 
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KNT42 
KNT45 
KNT47 


KNT36 
KNT40 


KNT28 
KNT31 
KNT35 


KNT23 
KNT25 


KNT20 
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Figure 1. Annual incidence of pertussis in Korea and list of tested strains in this study. Data of annual incidences in Korea were 
retrieved from the WHO data and disease web statistics system of KCDC. The upper line and number indicate annual pertussis 
incidence of Korea from 2000. No isolate tested in this study is shown in the lower table. Old isolates are defined as strains isolated 
before 1999, as indicated. The strain number of the isolates before 1980 was assigned as K00 and the strains named as KNT00 were 
the isolates after 1990. 



Bordetella MLST sequence-type database. For antigenic 
determinant genes the reported DNA sequences of each 
genotype were compared with those of the tested strains. 
Table 2 summarizes variation sites of 10 genes. Except 
Prn and TcfA the variation types of the other eight genes 
were point mutations at a single indicated site. However, 
Prn and TcfA genes showed systemic deletion of repeat 
units in the indicated region. To determine the genotype 
of each gene, reported DNA sequences were retrieved 
from the NCBI GenBank and aligned with tested 
sequences by the multiple alignment method using the 
MEGA program [26]. From this aligned DNA matrix the 
genotypes of tested strains were confirmed by compar- 
ison with the variable sites. If new variation sites were 
observed, the sequence was defined as new genotype 
and submitted to NCBI GenBank. 

For sequence types of housekeeping genes, genotype 
profiles of each isolate were queried to Bordetella 
MLST sequence-type database. For antigenic determi- 
nant genes we independently assigned STs according to 
the frequencies of genotype profiles confirmed in this 
study. 

The relatedness among each ST was analyzed by the 
START program package [27]. From this analysis, 
a clustering result was represented as a dendrogram by 
the Unweighted Pair Group Method with Arithmetic 
Mean (UPGMA) method. In addition, we compared the 
two dendrograms produced by analysis of housekeeping 
genes and antigenic determinant genes. The clonality of 



each ST was also analyzed by BURST algorithm and 
represented as a minimum spanning tree. 

3. Results 

3.1. Genotype polymorphisms in housekeeping 
genes and antigenic determinant genes 

Briefly, the polymorphic level of antigenic determi- 
nant genes was higher than housekeeping genes. As 
indicated in Table 3, the Korean genotypes of 7 
housekeeping genes showed similar patterns of reported 
global genotypes [28]. Except Adk and Tyr B genes, only 
single a genotype was observed in the other five genes. 
There were two different genotypes in Adk and Tyr B 
genes; the major type was type 1 and type 3, respec- 
tively. However, in the case of antigenic determinant 
genes a more complicated genotype pattern was 
observed. Currently, the reported global genotypes were 
45 in 10 tested genes. From these reported types 23 
genotypes (52.3%) were confirmed in the Korean 
isolates (Table 3). Under the consideration of average 
genotype per gene, antigenic determinant gene and 
housekeeping gene showed 2.2 and 1.3, respectively, 
suggesting that while the genetic variations were rare 
events in housekeeping genes, more genetic variations 
occurred in antigenic determinant genes. 

From the 10 antigenic determinant genes global 
genotype variations were mainly reported in PtxSl, Prn, 



00 

Table 1. PCR primers used in this study 

Gene Size of gene Primer sequence Size of analyzed region Reference 

Housekeeping genes 



Adenylate kinase (Adk) 


657 bp 


F : 
R 


AGCCGCCTTTCTCACCCAACACT 
TGGGCCCAGGACGAGTAGT 


530 


bp 


(30 bp- 


560 bp) 


Bordetella MLST sequence typ< 
database (www.Pubmlst.org) 


Fumarate hydratase class II (FumC) 


1392 bp 


F : 
R 


CGTGAACCGGGGCCAGTCGTC 
GGCCAGCCAGCGCACATCGTT 


498 


bp 


(396 bp- 


-894 bp) 




Serine hydroxymethyl-transferase (GlyA) 


1248 bp 


F : 
R 


CAACCAGGGCGTGTACATGGC 
CCGCGATGACGTGCATCAG 


499 


bp 


(402 bp- 


-901 bp) 




Aromatic amino-acid aminotransferase (TyrB) 


1203 bp 


F : 
R 


CGAGACCTACGCTTATTACGAT 
TGCCGGCCAGTTCATTTT 


487 


bp 


(462 bp- 


-949 bp) 




Isocitrate dehydrogenase (led) 


1257 bp 


F : 
R 


CTGGTCCACAAGGGCAACAT 
ACACCTGGGTGGCGCCTTC 


512 


bp 


(704 bp- 


-1216 bp) 




Cytosol aminopeptidase (PepA) 


1500 bp 


F : 
R 


CGCCCCAGGTTGAAGAAAATCGTC 
ATCAGGCCCACCACATCCAG 


508 


bp 


(454 bp- 


-962 bp) 




Phosphoglucomutase (Pgm) 


1383 bp 


F : 
R 


CGCCCATGTCACCAGCACCGA 
CGCCGTCTATCGTAACCAG 


548 


bp 


(672 bp- 


-1216 bp) 





Antigenic determinant genes 



Pertussis toxin subunit l (PtxSl) 


8l0bp 


F : 
R 


ACCATCAAAACGCAGAGGGGAAGA 
AGGAGGGCCAACGGCAGA 


935 bp (-29 bp to +98 bp) 


[25] 




Pertactin (Prn) 


2739 bp 


F : 
R 


GCCAATGTCACGGTCCAA 
CGGATTCAGGCGCAACTC 


1,410 bp (649 bp-2059 bp) 


[23] 




Filamentous hemagglutinin (FHA) 


10,772 bp 


F : 


TCGCCATTTCGGCGCACG 


372 bp (2348 bp-2720 bp) 


This study 






R 


AGATCGAGCTGCGCGCCG 








Fimbriae 2 (Fim2) 


624 bp 


F : 
R 


GCATAGACCATCTTGTATG 
CCGGCCGGGCTCCTTGAG 


388 bp (334 bp to +98 bp) 


[12] 




Fimbriae 3 (Fim3) 


615 bp 


F : 
R 


CCCCCGGACCTGATATTGTGATG 
GCTGAGCGTGCTGAAGGACAAGAT 


810 bp (-68 bp to +119 bp) 


[24] 




Outer-membrane protein Q (OmpQ) 


1095 bp 


F : 
R 


ATGCGTCGTCTTCTCGTC 
TCAGAAGCGCTGGGTCAT 


1,095 bp (1 bp- 1,095 bp) 


This 


study 


Bordetella antitransporter-protein C (BapC) 


2268 bp 


F : 
R 


GACAACGGTGTCTGGGGC 
GCGCAGGTGGAACGTCCA 


935 bp (1,327 bp-2,262 bp) 


This 


study 


Adenylate cyclase toxin (CyaA) 


5121 bp 


F : 
R 


ACCGCCTACGGCAAGCGC 
GCCGCCTTCAAGGGTATC 


1,043 bp (2,935 bp-3,987 bp) 


This 


study 


Virulence-activated gene (Vag8) 


2748 bp 


F : 
R 


ATGGCAGGACAAGCGAGG 
CCGCGTACCCGTCAACGT 


741 bp (1 bp-741 bp) 


This 


study 


Tracheal colonizing factor (TcfA) 


2019 bp 


F : 
R 


ATGCACATTTACGGAAATATGA 
TATGCGTGCCCGGGTCATAG 


1,105 bp (1 bp-1,105 bp) 


This 


study 



ZJ 
OQ 

fl) 
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Table 2. Genotype polymorphism of antigenic determinant genes analyzed in this study 



Gene 



Polymorphic sites 



Allele 


rvClCI C1JVCS> 


PtxSlA(l) 


AJ245366 


PtxSlB(2) 


AJ245367 


PtxSlC(3) 


[25] 


PtxSlD(4) 


AJ245368 


PtxSlE(5) 


AJ006151 


PtxSlF(6) 


AJ506994 


PtxSlG(7) 


a t r r\ s r\c\ c 

AJ506995 


PtxSlH(8) 


HM2 12424 


Prn(l) 


AF456359 


Prn(2) 


AF348484 


Prn(3) 


AF48485 


Prn(4) 


AJ011015 


Prn(5) 


AJ011016 


Prn(6) 


AF456357 


Prn(7) 


AJ133748 


Prn(8) 


AJ133245 


Prn(9) 


AF456356 


r Lily 


A T7R487S 


Prn(ll) 


AJ507642 


FHA(l) 


X52156 


FHA(2) 


AJ420989 


Fim2(l) 


AY845256 


Fim2(2) 


AJ420988 


Fim3(l) 


X51543 


Fim3(2) 


AY464179 


Fim3(3) 


AY464180 


Fim3(4) 


AY464181 


Fim3(5) 


AY845257 


OmpQ(l) 


U16266 


OmpQ(2) 


AJ420989 


BapC(l) 


AF081499 


BapC(2) 


[12] 


CyaA(l) 


Y00545 


CyaA(2) 


BX470248 



Pertussis toxin subunit 1 
(PtxSl) 



Pertactin (Prn) 



Filamentous 
hemagglutinin (FHA) 



Fimbriae2 (Fim2) 



Fimbriae3 (Fim3) 



Outer membrane protein 
Q (OmpQ) 



Bordetella autotransporter 
protein C (BapC) 

Adenylate cyclase toxin 
(CyaA) 



Tracheal colonizing factor 
(TcfA) 



585 600 681 691 



II 



.1 69 

IL 



810 

ZD 



C/A C/A C/T 



C/T C/T G/A G/A 



1 828 

DDD[ 

Systemic deletion of repeat unit 



2763 



:ddo 



]DDf 



10780 

]DDD 



260 263 320 388 



615 

ZD 



C/T 



C/A C/T G/A G/A 



1025 1 09 5 
I 



C/T 



DDDL 



2280 

1 



DDDC 



G/C A/G 



5121 

innn 



Systemic deletion of repeat unit 



□uuu 



Vag8(l) 


U90124 


Vag8(2) 


AJ420992 


TcfA(l) 


U16754 


TcfA(2) 


AJ009785 


TcfA(3) 


AJ420991 


TcfA(4) 


AJ507643 


TcfA(5) 


AJ420992 


TcfA(6) 


AY375533 


TcfA(7) 


AM238667 



120 



S.-O. Jung, et at 



Table 3. Summary of genotype and frequency observed in Korean isolates of B. pertussis 



Gene No. reported genotypes No. confirmed genotypes Type and frequency (%) 

Housekeeping genes 

Adk 2 2 Type 1 (96.7), type 2 (3.3) 

Fum C 1 1 Type 1 (100) 

Gly A 1 1 Type 1 (100) 

Tyr B 2 2 Type 1 (9.8), type 3 (90.2) 

led 1 1 Type 1 (100) 

Pep A 1 1 Type 1 (100) 

Pgm 1 1 Type 1 (100) 



Antigenic determinant genes 

PtxSl 8 5 Type 1 (93.2), type 2 (1.4), type 4 (1.4), 

type 5 (2.7%), type 8 (1.4) 



Prn 


11 


2 


Type 1 (96.7), type 6 (3.3) 


Fha 


2 


2 


Type 1 (96.7), type 2 (3.3) 


Fim2 


2 


1 


Type 1 (100) 


Fim3 


5 


2 


Type 1 (21.3), type 2 (78.7) 


Omp Q 


2 


2 


Type 1 (3.3), type 2 (96.7) 


Bap C 


2 


2 


Type 1 (96.7), type 2 (3.3) 


Cya A 


2 


2 


Type 1 (3.3), type 2 (96.7) 


Vag 8 


2 


2 


Type 1 (3.3), type 2 (96.7) 


Tcf A 


9 


2 


Type 1 (3.3), type 2 (96.7) 



Adk, adenylate kinase; fumC, fumarate hydratase class II, glyA, serine hydroxymethyltransferase; tyrB, aromatic amino-acid aminotransferase; icd, 
isocitrate dehydrogenase; pep A, cytosol aminopeptidase; pgm, phosphoglucomutase; ptxSl, pertussis toxin; prn, pertactin; fha, filamentous hemag- 
glutinin; fim2, fimbriae2; fim3, fimbriae3; ompQ, outer-membrane protein Q; bapC, Bordetella autotransporter protein C; cyaA, adenylate cyclase toxin; 
vag8, virulence-activated gene; tcfA, tracheal colonizing factor. 



Fim3, and TcfA genes (Table 3). However, in the Korean 
isolates, polymorphisms of Fim3 and TcfA were rela- 
tively low as compared with global genotypes. The most 
polymorphic gene was PtxSl; 5 types (1, 2, 4, 5, and 8) 
were observed and the major form was type 1 (93.2%). 
Especially, the PtxSl genotype of KNT 72 was 
confirmed as a new genotype. This variation was a silent 
point mutation occurring at site 600 bp where base C 
was changed to T. We listed up this new sequence to 
GenBank (HM2 12424) and assigned this genotype as 
type 8. 

Compared of the prevalence of PtxSl genotype from 
previous reports of other countries [10,12,16,29—33] 
confirmed that the most frequent global genotype was 
type 1, except in Japan (Table 4) and the next most 
frequent form was type 2. Other genotypes (types 3—8) 
were rarely reported. Especially, in Australia and Italy 
only type 1 was reported. However, Japan and the USA 
showed different patterns; in the former, type 2 showed 
higher frequency than type 1, whereas in the latter the 
ratio of type 2 was higher than in other reported coun- 
tries. In the case of the Korean isolates the major 
genotype was type 1 similar to other countries, but rare 
types (type 4 and 5) were also observed, and even a new 
type (type 8). Therefore Korean isolates showed the 
most diverse pattern in the PtxSl gene among countries. 

While the PtxSl genotype was most diverse in the 
Korean isolates, the Prn genotype showed a different 
pattern from those in the reports by other countries 
[10,12,16,29—33]. Currently, 11 different prn genotypes 



were reported and the major difference between geno- 
types was attributed to deletion of repeated amino acid 
sequence units (GGxxP and PQP) [13]. From these 
genotypes the frequently observed genotypes were type 
1, 2, and 3. Other genotypes (types 4—11) were rarely 
observed. As indicated in Table 5, more than types were 
reported from other countries, and the major genotype 
was also different in each country. Japan, the UK, and 
Australia showed type 1 as the major genotype, whereas 
type 2 was major genotype in Taiwan, Finland, and 
Sweden. In the Korean isolates, 3 genotypes were 
observed; the major genotype was type 1. Notably, the 
rare type 6 was also observed. Therefore, the Prn gene 
of Korean isolates was not as diverse a gene as those 
reported from other countries. 

3.2. Multilocus sequence type determination 
and comparative analysis of sequence types 

Currently, a total of 43 STs (housekeeping gene 
sequence types; HSTs) for seven housekeeping genes 
were reported in Bordetella MLST database (http:// 
pubmlst.org). These HSTs were established from anal- 
ysis of collected genotype data from three Bordetella 
species (B. pertussis, Bordetella bronchiseptica, and 
Bordetella parapertussis) by Diavatopoulos et al. [28]. 
From these HSTs, only three HST types (1,2 and 24) are 
originated from B. pertussis species. According to the 
analysis of clonality by Bordetella MLST database, 43 
STs are reorganized into four clonal complexes. From 



Table 4. Comparison of PtxSl genotypes with the reports from other countries 



Allele type 


GenBank 
accession 
no. or ref. 


Korea 
(n = 81 b ) 1} 
lyoo — zuiu 


Japan 

(n = 107) 2) 

iqoq onm 
l^oo — ZUU1 


Taiwan 
(n = 80) 3) 

1 QQO 1(\(\A 

lyyo — ZUU4 


USA 
(n = 152) 4) 

inn/: 


Finland 

(n = 122) 5) 


UK 

(n = 335) 6) 

lyZV — ZUUZ 


Sweden 
{n = 1006) 7) 

ly IK) — ZUUj 


Australia 
(n = 46) 8) 

ly IK) — ZUUj 


Italy 

(n = 30) 9) 

1 QQ^ 1 QQ^ 

lyyj — lyyj 


S1A (1) 


AJ245366 




23 


79 


98 


115 


288 


998 


46 


30 


SIB (2) 


AJ245367 


6 


84 


1 


51 


7 


36 


8 






SIC (3) 


[25] 




















SID (4) 


AJ245368 


1 






3 












S1E (5) 


AJ006151 


2 


















S1F (6) 


AJ506994 




















S1G (7) 


AJ506995 




















S1G (8) 


HM2 12424 


1 


















"Strain number reconstructed from data of the indicated reference; b Number of tested strain was reconstructed from cited 
Cassisay et al, 2000 [10]; 5) Elomaa et al, 2005 [31]; 6) Packard et al, 2004 [12]; 7) Hallander et at, 2005 [16]; 6 


references (superior number). 1) In this study; 2) Kodama et al, 2004 [29]; 3) Yao et al, 2005 [30]; 4) 
) Byrne and Slack, 2006 [32]; 9) Mastrantonio et al, 1999 [33]. 


Table 5. Comparison of Prn genotypes with reports from other countries 














Allele type 


GenBank 
accession 
no. or ref. 


Korea 
(n = 81 b ) 1} 
1968-2010 


Japan 
(n = 107) 2) 
1988-2001 


Taiwan 
{n = 80) 3) 
1998-2004 


USA 
(n = 152) 4) 
1935-1996 


Finland 
(n = 122) 5) 
1953-2003 


UK 
{n = 335) 6) 
1920-2002 


Sweden 
(n = 919) 7) 
1992-2000 


Australia 
{n = 46) 8) 
1970-2003 


Italy 
(n = 129) 9) 
1993-1995 


Prn (1) 


AF456359 


72 a 


88 


6 


86 


26 


213 


51 


38 


8 


Prn (2) 


AF348484 


7 


18 


72 


66 


92 


111 


748 


3 


53 


Prn (3) 


AF.48485 




1 


2 




2 


6 


116 


5 


65 


Prn (4) 


AJ011015 










2 










Prn (5) 


AJ011016 


















3 


Prn (6) 


AF456357 


2 


















Prn (7) 


AJ133784 




















Prn (8) 


AJ1 33245 




















Prn (9) 


AF456356 




















Prn (10) 


AJ784875 




















Prn (11) 


AJ507642 





















a Strain number reconstructed from data of the indicated reference; Number of tested strain reconstructed from cited references (superior number). 1) In this study; 2) Kodama et al, 2004 [29]; 3) Yao et al, 2005 [30]; 4) 
Cassisay et al, 2000 [10]; 5) Elomaa et al, 2005 [31]; 6) Packard et al, 2004 [12]; 7) Hallander et al, 2005 [16]; 8) Byrne and Slack, 2006 [32]; 9) Mastrantonio et al, 1999 [33]. 
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Table 6. Summary of genotype profiles from Korean isolates of B. pertussis 



Genotype 


ST of Housekeeping genes (HST) GlyA 


Adk 




led 


Tyr B 


Pep A 


Pgm 




Fum C 


1 (4.9%) 1 


1 




1 


1 


1 


1 




1 


2 (92.6%) 1 


1 




1 


3 


1 


1 




1 


24 (2.5%) 1 


2 




1 


1 


1 


1 




1 


Genotype 


ST of Antigenic determinant genes (AST) PtxSl 


Prn 


FHA 


Fim2 


Fim3 OmpQ 


BapC 


CyaA 


Vag8 


TcfA 


1 (72.8%) 1 


1 


1 




2 2 


1 


2 


2 


2 


2 (8.6%) 1 


2 


1 




1 2 


1 


2 


2 


2 


3 (6.2%) 1 


1 


1 




1 2 


1 


2 


2 


2 


4 (2.5%) 5 


6 


2 




1 1 


2 


1 


1 


1 


5 (1.2%) 4 


1 


1 




2 2 


1 


2 


2 


2 


6 (7.4%) 2 


1 


1 




1 2 


1 


2 


2 


2 


7 (1.2%) 8 


1 


1 




2 2 


1 


2 


2 


2 



Adk, adenylate kinase; fumC, fumarate hydratase class II, glyA, serine hydroxymethyltransferase; tyrB, aromatic amino-acid aminotransferase; icd, 
isocitrate dehydrogenase; pep A, cytosol aminopeptidase; pgm, phosphoglucomutase; ptxSl, pertussis toxin; prn, pertactin; fha, filamentous hemag- 
glutinin; fim2, fimbriae2; fim3, fimbriae3; ompQ, outer-membrane protein Q; bapC, Bordetella autotransporter protein C; cyaA, adenylate cyclase toxin; 
vag8, virulence-activated gene; tcfA, tracheal colonizing factor. 



Table 7. Lineage assignment analysis by genotype profiles of housekeeping genes and antigenic determinant genes 



Gene type 



Group a 



ST Frequency SLV DLV SAT 



Housekeeping genes 




l b 


4 


2 




2 


75 


1 


1 


24 


2 


1 


1 



Antigenic determinant genes 




6 


6 


1 


4 


0 


l b 


59 


3 


2 


0 


3 


5 


3 


2 


0 


5 


1 


2 


2 


1 


7 


1 


2 


2 


1 


2 


7 


1 


2 


2 


4 


2 









Singleton 



a Group definition (N—l); b Central types. SLV, single locus variation; DLV, double locus variation; SAT, satellite strain. 
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these clonal complexes, clonal complex II is the group 
of B. pertussis comprising HSTl, HST2 and HST24. 

In the Korean isolates (Table 6), 3 ST types (HSTl, 
HST2, and HST24) were also observed and the 
frequencies of each HST were 92.6% for HST2, 4.9% 
for HSTl, and 2.5% for HST24 as similar to the pattern 
of MLST database. Therefore the distributing strains in 
Korea were in a single clonal complex II, as indicated in 
the MLST database. 

In the case of STs of antigenic determinant genes 
(antigenic determinant gene sequence types; ASTs), 
a more divergent pattern was observed. The Korean 
isolates were divided into seven ASTs based on the 
genotype profiles of 10 antigenic determinant genes. 
Because there was no database for ASTs of 10 antigenic 
determinant genes, we independently assigned ASTs 
according to the observed frequencies. As shown in 
Table 6, type 1 (most frequent form in the Korean 
isolates) showed 72.8% of frequency and other types 
(2—7) showed frequencies <10%. 

When we analyzed relatedness among these HSTs by 
BURST algorithm using the START program package 
[27], the STs confirmed in the Korean isolates were 
grouped in single clonal complex as expected and the 
central ST type was determined as HSTl (Table 7). The 
strains included in HSTl were Kl, K2, K5, and K23. 
These isolates were old strains isolated in 1968 and 1972 
(Figure 1). In the old strains, the frequently observed ST 



profile was HSTl, but the major sequence type was 
changed to HST2 in the strains isolated from 2000. 
Although HST24 was also observed, the frequency was 
too low and only two strains were transiently observed 
in 2000 and 2001. Therefore a major sequence-type 
transition from HSTl to HST2 was confirmed accord- 
ing to time period in the Korean isolates. 

In the AST analysis, as shown in Table 7, it was 
revealed that the Korean isolates were categorized into 
two clonal complexes, except AST4. AST4 showed 
a unique and quite different genotype as compared with 
other genotype profiles. This sequence type was not 
a frequently observed form in the Korean isolates, and 
only two strains (KNT22 and KNT23) were observed in 
2000—01. Therefore it was inferred that this sequence 
type was transiently introduced type from outside Korea. 

These two clonal complexes were assigned as anti- 
genic clonal complex (ACC) I and II. ACC I was 
composed of AST2, AST3, and AST6. The property of 
this clonal complex was old or new sequence type. 
AST3 was the central type of ACC I and emerged in 
2000 and reemerged in 2009. AST6 was only confirmed 
in old strains isolated in 1968 and 1972. However, AST2 
was the new sequence type that appeared in 2009, and 
the major difference was a change of genotype of Prn 
from 1 to 2. 

The other clonal complex, ACC II, was composed of 
AST1, AST5, and AST7. The major sequence type of 



Housekeeping gene 



Antigenic determinant gene 



HSTl 



:K1, K2, K5, K23. 



HST2 



K3, KNT 40 



:KNT 75 , 76 . 78 , 80, 82 , 83 , 84 



KNT 8, 16, 17, 73, 81 



Other 59 strains 



KNT 72 



KNT 39 



:AST6 



*AST2 



:AST3 



:asti 



HST24 



KNT 23. KNT 22 



IAST7 
•AST5 

1 AST4 



I 1 

0.02 



0.1 



Figure 2. Comparison of UPGMA trees originating from analysis of housekeeping genes and antigenic determinant genes. Two 
UPGMA trees of genotype profiles deriving from housekeeping genes (HSTs) and antigenic determinant genes (ASTs) represented 
as merged form. The left-side tree was derived from housekeeping genes and the right side tree from antigenic determinant genes. 
The line box indicates sequence type clustered by housekeeping genes analysis and dot box strain sequence type clustered by 
antigenic determinant genes analysis. 
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Table 8. Transition pattern according to time period of major STs in housekeeping genes and antigenic determinant genes 



Time periods 




Housekeeping genes 




Antigenic determinant 


genes 


Major ST 


Genotype profile a 


Frequency (%) 


Major ST 


Genotype profile b 


Frequency (%) 


Tohama I 


1 


1,1,1,1,1,1,1 






2,1,2,1,1,2,1,2,2,2 




Old (1968/1975) 


1 


1,1,1,1,1,1,1 


80 


6 


2,1,1,1,1,2,1,2,2,2 


100 


1999 - 2008 


2 


1,1,1,3,1,1,1 


96.9 


1 


1,1,1,1,2,2,1,2,2,2 


87.5 


Outbreak in 2009 


2 


1,1,1,3,1,1,1 


100 


2 


1,2,1,1,1,2,1,2,2,2 


58.3 



a Gene order of housekeeping genes: GlyA, Adk, led, TyrB, PepA, Pgm, FumC; b Gene order of antigenic determinant genes: PtxSl, Prn, FHA, Fim2, Fim3, 
OmpQ, BapC CyaA, Vag8, TcfA. 



ACC II was AST1 and AST1 was also determined to be 
the central type. The strains of AST1 were frequently 
observed from 1999 to 2008. However, the frequency 
was greatly decreased in 2009. 

These different results between HSTs and ASTs 
indicate that there were different driving forces for 
genetic events. This was clearly confirmed when these 
two gene categories were compared by UPGMA tree. As 
shown in Figure 2, all the strains belonging to HST2 
showed more diverse genotype variations in antigenic 
determinant genes, and finally diverged into six ASTs, 
indicating that there were different evolutional changes, 
such as vaccination in antigenic determinant genes 
[9,11,25,30,33-36]. 



4. Discussion 

The emergence mechanisms of variants are not 
exactly known. However, it is considered that there are 
two major ways for variants to emerge. One is that 
existing strains acquire different genotype and pheno- 
type by mutation; the other is that distributing strains are 
selected from the population by factors such as vacci- 
nation. From these two considerations, it is not clear 
which is predominant. In the case of B. pertussis, this 
species is known to be a homologous population by 
MLEE pattern analysis [37]. The variant sequence was 
also rarely confirmed in housekeeping genes [11]. 
However, there were many reports about genotype 
changes in antigenic determinant genes. Especially it 
was also suggested by Smith et al. [4] that because B. 
pertussis is an exquisitely adapted pathogen to the 
human respiratory tract by evolutionary process, its 
rapid evolution rate should be considered in developing 
a vaccine and protecting efficacy. Therefore it is inferred 
that emergence of the genotype variants in antigenic 
determinant genes is the active process of B. pertussis to 
avoid the protective effect of vaccination. 

These genotype variations also showed a transition 
pattern according to the vaccine introduction. In the 
Netherlands, the UK, Sweden, Taiwan, and Japan 
[9,1 1 — 14,25,32] it was reported that the major sequence 
types before vaccine introduction gradually changed 
thereafter. That is, the distributed strains containing the 



same genotype with the vaccine strain were gradually 
reduced after vaccination, while other strains containing 
different genotypes appeared. In Korea, similar geno- 
type changes were confirmed by MLST analysis. One of 
the advantages of MLST is that it is easy to see transi- 
tional trends of bacterial type over the long term. In this 
analysis, the Korean isolates also showed a transitional 
pattern in antigenic determinant genes. The first transi- 
tion was observed in the isolates of 1999. In this stage, 
the genotype profiles of the major isolates were changed 
from AST6 (2,1,1,1,1,2,1,2,2,2) to AST1 
(1,1,1,1,2,2,1,2,2,2). As indicated in Table 8, this was 
attributed to genotype change of PtxSl and Fim3 
genes. In addition, major ST of housekeeping gene also 
changed from HST1 (1,1,1,1,1,1,1) to HST2 
(1,1,1,3,1,1,1). When we considered that the acellular 
vaccine derived from Tohama I strain was introduced 
from 1982 [19], these genotype change patterns also 
showed a tendency to avoid vaccine- derived genotypes 
as described above. 

A second transition was also observed, in 2009. At 
that time, pertussis incidence was greatly increased by 



ACCI I 



.Qld strain 
AST6 




Figure 3. Minimum spanning tree of Korean isolates by 
MLST profiles of antigenic determinant genes. A minimum 
spanning tree was constructed by sequence type profiles 
deriving from antigenic determinant genes. Lines-circles 
indicate sequence type; numbers in parentheses are the number 
of strains included in each sequence type. Dots-circles indicate 
the clonal complex confirmed by BURST analysis. Because 
the distance between AST4 with AST3 was too long, we 
omitted (asterisk). AST2 was newly appearing sequence type 
when the pertussis incidence was greatly increased. 
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500% as compared with the average for the previous 
decade. AST of the isolates in 2009 changed from AST1 
to AST2 (1,2,1,1,1,2,1,2,2,2). This ST change was 
attributed to genotype changes of Prn and Fim3 genes. 
Especially, the genotype change of Prn gene was first 
confirmed in the Korean isolates and new ST was 
diverged. It was also clearly confirmed by minimum 
spanning tree (Figure 3). As shown in Figure 3, the 
ASTs of the Korean isolates were divided into two 
clonal complexes, ACC I and ACC II. From these clonal 
complexes it was revealed that AST1 of ACC II was 
distributed until 2008 and AST2 of ACC I newly 
appeared from 2009. 

Genotype changes of the distributing strains are still 
actively progressing events. Although the polymorphic 
status of Korean isolates was relatively lower than other 
countries, there is enough evidence to warn of the risk of 
pertussis outbreak in Korea. For this reason, construc- 
tion of an efficient and continuous surveillance system 
to detect emergence of genotype variants is required. 
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